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We report on the status of an ALPHA Collaboration project to extract quantities for B physics 
phenomenology from Nf = 2 lattice simulations. The framework is Heavy Quark Effective The- 
PsJ ory (HQET) expanded up to the first order of the inverse b-quaik mass. The couplings of the 

effective theory are determined by imposing matching conditions of observables computed in 
HQET with their counterpart computed in QCD. That program, based on Nf = 2 simulations in a 
small physical volume with Schrodinger functional boundary conditions, is now almost finished. 
On the other side the analysis of configurations selected from the CLS ensembles, in order to 



5—1 measure HQET hadronic matrix elements, has just started recently so that only results obtained 

at a single lattice spacing, a = 0.07 fm, will be discussed. We give our first results for the b-quaik 
mass and for the B meson decay constant. 
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1. Introduction 

Now that the amount of available experimental data on beauty physics coming from the B 
factories and the Tevatron is large, and even more data are expected from the LHC, precision 
tests of the Standard Model and searches for New Physics have become possible in this area of 
flavor physics. Unfortunately, the theoretical uncertainty, mainly because of the difficult to estimate 
long-distance effects due to confinement, is currently limiting the impact of future experimental 
measurements on New Physics models. Lattice QCD makes it possible to reach a few percents 
of theoretical error on those non-perturbative hadronic contributions, but care is needed to obtain 
reliable results for ft-quark physics. Indeed one has to keep under control simultaneously the finite 
size effects and, particularly, the discretisation effects (the lattice spacing should be smaller than 
the Compton wavelength of the Z?-quark) induced by the simulation. In practice it is not possible 
to control both effects in one simulation. Different approaches have been proposed in the literature 
(see for example [1] for a recent review). The ALPHA Collaboration has followed a strategy 
discussed in detail in [2] - [4]: it is based on the use of HQET, in which the hard degrees of 
freedom ~ mb are integrated out and taken into account by an expansion in the inverse Z^-quark 
mass /Mb . As discussed in those papers and also in earlier work, the benefit is the suppression of 
large discretisation effects which may arise in hadronic quantities when the theory is regularised on 
the lattice. The difficult aspect of that method is that a matching with QCD, which is the field theory 
believed to describe the strong interactions, is needed to absorb ultraviolet divergences appearing 
in the effective theory. In HQET those come as inverse powers of the lattice spacing and thus have 
to be removed non-perturbatively before the continuum limit can be taken. This method has been 
tested in the quenched approximation to study the meson spectrum [5] and to determine the 
ft-quark mass [3], the decay constant /b, [6], and the coupling gB*B7t [7]. 

The strategy, sketched in Fig. 1 , is now being extended to the more realistic Nf = 2 situation. 
In the next section, more details are given about the determination of HQET couplings, performed 
in a small physical volume (of space extent ~ 0.5 fm). In section 3, we report on the analysis of a 
subset of ensembles generated within the CLS effort, to get HQET energies and matrix elements. 
Section 4 contains our conclusions. 

2. Computation of the relevant HQET parameters 

We write the HQET Lagrangian as 



^HQEt(-^) = ^stat(-^) - Oh:m^\;m{x) - (0^pm^spm{x) , 



(2.1) 



where the lowest order (static) term is 



(2.2) 



and the first order corrections in 1 / m\j are 



(2.3) 
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We are also interested in the time component of the heavy-light axial current Aq. Considering only 
the terms which contribute to zero-momentum correlation functions, we write it as 



.HQET 



(x) 



(2.4) 
(2.5) 



^HQET^^O) ^j^j nibare (the energy shift 



where Vf denotes the symmetric derivative. 

In this section we present the computation of COtin, G)spin,^A 
which in the static theory absorbs the I /a divergence of the static energy and at order 1 / nib absorbs 
a 1/a^ term). These parameters can be used for a computation of the Z^-quark mass, the heavy-light 
meson decay constants (/b or /bJ, as well as for a determination of the spectrum of heavy-light 
mesons, including the hyperfine mass splitting. We follow the strategy presented in [4], where the 
computation was done in the quenched approximation. We recall here the essential ingredients and 
refer the reader to this paper for any unexplained notation and for more detailed explanations. This 
computation is done non-perturbatively (in the strong coupling), and at the 1/nib order of the heavy 
quark expansion. The light quarks are simulated with 2-flavor Clover-improved Wilson fermions, 
and for the discretisation of the heavy quark we use the so-called HYPl and HYP2 actions [8]. 
More details about the implementation can be found in [9]. The simulations considered in this 
section use the Schrodinger functional setup. 

Following [4], we define a set of observables <I>,=i .5 that we collect in a vector <I>. In the 
continuum and large volume limits, is proportional to the meson mass and <I>2 to the logarithm 
of the decay constant, respectively. is used to determine the counter-term of the axial current, 
and <I>4 5 for the determination of the kinetic and magnetic term, respectively. We first consider 



QCD 



HQET 




Si 



Figure 1: Sketch of the strategy followed by the ALPHA Collaboration to compute B physics observables 
on the lattice. 
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a small volume (of linear space extent Li ~ 0.5 fm), where the ft-quark mass can be simulated 
with discretisation effects under control. In this volume, we compute the five observables <J>,- in 
QCD, for four different lattice spacings. In this set of simulations, the light quark masses are 
set to 0, while for the (RGI) heavy quark mass M we have chosen nine different values such 
that z = L\M = 4,6,7,9, 11, 13, 15, 18,21. The lightest masses correspond approximately to the 
charm and the heaviest to the bottom quark. In this work we will focus on the heaviest ones. The 
continuum limit of each observable is obtained by a Unear extrapolation in {a/L\Y of the three 
finest lattice spacings (except for the two heaviest masses where only the two finest lattices are 
used): <I>|^*^'^(Li,M,0) = \ima^Q^^^{L\,M,a). This set of simulations is represented by in 
Fig. 1. The continuum extrapolations of the first two observables are shown in Fig. 2. 




(a/Li)2 xio" (a/Li)2 xio"^ 

Figure 2: Continuum extrapolation of the QCD observables <I>i and 4>2, which are proportional to the finite 
volume meson mass and to the finite volume logarithm of the decay constant, respectively. For <I>i we have 
included an error in the continuum coming from the renormalisation of the quark mass (cross on the left). 
Results are shown for the nine different values of the heavy quark mass. 



In another set of simulations called 52, we compute the corresponding quantities in the effec- 
tive theory, using the same value of the physical volume. We then impose the QCD observables 
to be equal to their HQET expansion at the 1 / wib order ' . This matching can be written in the 
following way: 

a>Q™(Li,M,0) = T],(Z.i,fl) +£9,7(^1, a) (j)y(M,fl) , (2.6) 

where 77 and (p we. computed by lattice simulations for different values of the lattice spacing a. 
In other words the matching equations determine the set of pai^ameters w = (p^' [<I)*3'-'-' — r\]. For 
example, in the the static approximation (p is diagonal and (up to a factor L{) r\\ is given by the 
static energy, (p\ 1 is one and a)\ is the bare quark mass in static approximation. In table 1 , we list 
the various parameters together with their values in the static approximation and in the classical 
limit. 

We then compute T] and 9 in a larger volume of space extent L2 = 2L\ using the same set of 
lattice spacings as the one used in the previous step. This step is represented by the set ^3. The 
observables in this volume are obtained from the parameters «(M,a) determined in the previous 

'Terms of order l/"2^ are dropped without notice. 
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w, 


UClllllLlUll 


LidSSlCdl Value 


SLdLlC VdlUC 


ecu 

0)3 
0)4 
G)5 


'^^bare 
(Okm 


nib 


-l/(2nib) 
l/(2nib) 
l/(2nib) 


"Ce 

log(zr) 






Table 1: Notation for HQET parameters, their values in classical and static approximation. For the numerical value of 
c^"', we use the formula given in [10] . 

Step: 

^{L2,M,0) = lim{ri{L2,a) + q){L2,a)o){M,a)) . (2.7) 

The continuum limit can be taken because the divergences cancel out in the previous equation. This 
procedure can then be re-iterated until the volume reached is large enough for finite size effects to 
be negligible, typically around (2fm)^. In practice, it turns out that three different volumes are 
enough {L\,L2 and the large volume one). Thus, the HQET parameters that can be used in large 
volume simulations (denoted as ^4) are given by 

CO (M, a) = 9- i(L2, a) [<I>(L2,M,0)- 77 (L2, a)] . (2.8) 

Finally we perform a small interpolation in j8 in order to obtain the parameters at the lattice spacing 
used in the CLS ensembles. 

In general the equations (2.6) and (2.7) are meant to be taken at the 1/nib order of HQET, 
but they are of course valid in the static approximation if one sets the various pieces to their static 
values. In that case, the observables will be denoted by O**"" , and otherwise we write O^Q^^ 
also define the l/nib correction to be <I>(^/™) = O^Q^t _ <j)Stat illustration, we show some 

continuum extrapolations in the volume L2: the first two observables in the static approximation in 
Fig. 3, and the 1/nib corrections (^^^^"^^ and Oj'^™' (which are sensitive to the kinetic and magnetic 
parameter cOkin and ©spin) in Fig- 4. Note that the 1/nib terms are extrapolated linearly in the lattice 
spacing, whereas in the static case the continuum extrapolation is done quadratically in the lattice 
spacing. This is justified because we use an C?(a)-improved action and an C? (a) -improved static 
axial current. We end this section by collecting in Table 2 the values of the HQET parameters used 
in the numerical applications reported in the following. 

3. Extraction of HQET hadronic matrix elements 

For the computation of the large-volume hadronic matrix elements, we ought to take the con- 
tinuum limit and extrapolate in the up/down quark mass to the physical point where nij^ has its 
physical value. The continuum extrapolation will be left for future work, we here just take a quite 
small lattice spacing, a = 0.07 fm. We note that in our quenched computations the difference 
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2 4 6 '024 



(a/L2)2 xio" (a/L2)2 xio" 

Figure 3: Continuum extrapolation of the static approximation of <I>i and <I>2 in the volume of space extent 
L2- The two different colors correspond to two different discretisations of the static propagator : red for 
HYPl and blue for HYP2. Only the three finest lattice have been used in the continuum extrapolation. 




0.02 0.04 0.06 0.08 0.02 0.04 0.06 0.08 

a/L2 a/L2 



Figure 4: Same as Fig. 3 but for the l/mj, terms '™ and , which are chosen to be proportional to 

cOkin and (Ospm, respectively. 



between the continuum limit and a = 0.07 fm would be unnoticeable given the present errors. Con- 
cerning the extrapolation in the light quark mass, we take into account terms at NLO in the chiral 
expansion for the static approximation, but neglect terms of order ninght/nib- 
At the first order of the l/nn, expansion our main observables are given by 

'MB = mbare + E'''' + (0^^ E^"" + Wspi„ £^P^" , (3.1) 

mB-mB* = ^Wspin^'P'", (3.2) 

\og{a"^h./m^) = log(Z^Q^'^)+log(a^/2p^'^') + ^r«'"q 

-fcOkin/.'-" + «,pi„p^P'" + ci^'p^'" , (3.3) 

where b^^^ is an improvement coefficient (in practice we follow [10] for its numerical implemen- 
tation). The HQET energies and matrix elements have been measured on a subset of configuration 
ensembles produced within the CLS effort [11] with Nf = 2 flavors of C?(a)-improved Wilson- 
Clover fermions. We have collected in Table 3 the main characteristics of that subset. 
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i8 


LMq 




log(zr^) 


[log(Z^)](Vm) 


4'V« 




Wspin/a 


5.3 


15 


41.125(22) 


-0.120(28) 


0.039(26) 


-0.446(75) 


0.476(06) 


0.837(35) 



Table 2: HQET parameters determined with the standard set of matching conditions, that are used to obtain 
the central values of m\, and /b; the static quark action is HYP2 . 



3.1 Large volume techniques 

As interpolating fields we use quark bilinears 

Ok (x) = Wh (x) YqYs W*"^ (x) , 01 {x) = ^Iff'^ (x) 7o75 {x) , 



(3.4) 



built from the static quark field \l/h{x) and different levels of Gaussian smearing [12] for the light 
quark field 

\lfl''\x) = (1 + Koa^Af' Yi{x) , (3.5) 

with APE smeared links [13, 14] in the lattice Laplacian A, and with the same parameters as in [5]. 
We compute the following correlators 



Cf\t)=l^{0,ixo + t,y)0*{x))^^^^, 

^kin/spin^^^ = I(O,-(^0 + ?,y)O}(x)Oki„/sp:n(z] 



= lK^(xo+?,y)o*(x) 



(3.6) 
(3.7) 
(3.8) 



using stochastic all-to-all propagators to reduce the variance. From the N xN matrices of correla- 
tors C'"', c'''"/'P'", and C^\f), we solve the generalised eigenvalue problem (GEVP) 



C{t)Vn{t,to) = X„{t,to)C{tQ)Vn{t,to) , 

to get the energies and operators having the largest overlap with the state: 



aEf{tfy) = -In 



?in{t,tQ) 

0'{t)vi{t,to) 



\/<{t,to)Cij{t)vl{t,to) 



A„(fo + 2a,fo) 



t/2a 



(3.9) 

(3.10) 
(3.11) 



CLS 





a (fm) 


L^xT 


(MeV) 


# 


traj. sep. 


5.3 


0.07 


32^ X 64 


550 


152 


32 






32^ X 64 


400 


600 


32 






48^ x96 


300 


192 


16 






48^ x96 


250 


350 


16 



Table 3: Characteristics of the large volume simulations used so far to extract HQET energies and matrix 
elements. The last column is meant in terms of trajectories of length T = 1. 
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For a computation of the 1/mb corrections to static energies and matrix elements one needs to solve 
the GEVP at the static order only [15]. For instance, in the case of energies, one has 



Ef{t,to) = Ef''''\t,to) + (oET^"''\t,t^) + 0{(o') 



7eff,stat 



eff,l/m. 



aE, 



eff.stat 



In 



XT\t + aM) 



'J 



(3.12) 
(3.13) 

(3.14) 
(3.15) 



As illustrated in the left panel of Fig. 5, a comparison with quenched data indicates that the statis- 
tical eiTor on static energies is not an issue, even at small quark masses. Some cai^e has to be taken 
to control the contribution of excited states to extracted energies. For the energy levels determined 
from the GEVP, the leading corrections ai^e given by 



7eff,stat 



Stat -A£;v+l,/i? 



+ ... 



(3.16) 



where AE^,,! = E,n — En and the condition tQ>t/2 is necessary to prove eq. (3.16). In our analysis 
we have chosen a time range to extract the plateaux such that the corrections to Ef^^ are small 
compared to its statistical error; we found this to be the case for ?o > 0.3 fm. 

3.2 Zj-quarkmass 

Once the energies are obtained at the simulated dynamical quark masses, we still need to 
extrapolate them to the physical point. We have chosen pion masses below 450 MeV for the 



1.04 



0.96 



0.92 



--T-r^-.-f 



r f 



Nj = 2, m^ = 4O0 MeV 
Nj = 2, (11^ = 300 MeV 
N^ = 2, m^ = 25UMcV 
N =U, m =111 



0.4 



0.6 



0.8 
t[fm] 



1.2 



13.4 



2 13.2 



12.8 




fro™,)" 



Figure 5: Plateaux of static energies with A^=5 and fo = 5fl (left panel); the quenched data correspond to a 
lattice spacing quite close to the CLS set up at /3 = 5.3. We remind the reader that the absolute value of E**'^' 
has no real meaning without the subtraction of a linear I /a divergence. Right panel: chiral extrapolation of 
rQffig^' = '"o(£'*'^' +'«baie)- A large part of the errors (dashed part), originating from and Wbarej is common 
to all data points and therefore irrelevant for the quark mass dependence. 
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t [fm] (r„m^)' 

Figure 6: On the left panel we show the plateaux of cOspinB'*^" with A'=5 and fo — the quenched data 
correspond to a lattice spacing quite close to the CLS set up at /3 = 5.3. On the right panel we show the 
chiral extrapolation of ro{mQ — m^* ) with a LO fit. Here the dashed error due to COspin and ro is independent 
of the quark mass. 



extrapolation of all quantities reported in these proceedings. For the extrapolation of m^^' = + 
"^bare '^^^^ uscd the foUowing form 

romg^' = b + cml + dml , (3.17) 

where we once set d = 0, and once use d = —3rog^/{l6nf^) computed with the experimental 
value of and the recent lattice determination of ^ = 0.51(2) [7]. We show these fits in the right 
panel of Fig. 5. As our central value of romg'" at the physical point we take the average of the two 
extrapolations. Half the difference is included in the systematic errors. At present a much larger 
source of uncertainty is ro/a, whose determination is discussed in [17]. Moreover, an uncertainty 
of 5% has been added on the scale ro itself. We take ro = 0.475 ±0.025 fm [20-22]. Identifying the 
static approximation result, m^^^, with m^'^, we can determine the RGI Z?-quark mass from static 
HQET. Translating as usual to the MS scheme with 4-loop perturbation theory and the known 
A-parameter [18], we find 

m^{mf^y^%2 = 4.255(25),„(50)stat+renorm(?)« GcV , (3.18) 

where the first error comes from the uncertainty on ro, while the second error includes the statis- 
tical error on aE^^'^\ the uncertainty on the chiral extrapolation and the error on the quark mass 
renormalisation constant Zm in QCD [19]. The latter is currently dominating. The "(?)«" indicates 
that a continuum limit is not yet performed, but as mentioned earlier we expect only a smaller error 
due to that [4]. 

Amongst the 0{l/m\j) corrections to energies the hyperfine mass splitting is particularly in- 
teresting since here the 1 / nib term is the dominant one and is given by a single HQET parameter, 
OJspin- This contribution can hence be discussed separately. Its determination is encouraging, as 
far as the statistical uncertainty and the chiral extrapolation are concerned, as shown in Fig. 6. 
Plateaux have the same quality as in the quenched case and the statistical precision is good enough 
to perform a reasonable chiral extrapolation. Since we have seen in the quenched approximation 
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that cut-off effects can be sizable for this quantity, we do not quote the hyperfine splitting in MeV 
at the moment. 

Including the magnetic as well as the kinetic correction, we compute mg^^^ = wig^'^^''™, per- 
form the chiral extrapolation as for m^^\ and from equating mg^^^ = m^'^ we obtain the Z?-quark 
mass up to tiny A^/m^ effects as 



HQET 

^^^^ =4.276(25),o(50)stat+re„orm(?)« GcV . (3.19) 



Sources of errors are the same as in the static case. Some previous determinations read 

HQET 

^MS^MS) = 4. 320 (40), (48) GeV [3], 

Nf— 

^MS(,„MS) ^ 4.163(16) GeV [23]. 

In the future, a nice check of our result will consist in doing the same computation in a partially 
quenched set up; there, the experimental inputs will be the spectrum, and the determination of 
the hopping parameter of the strange quark will be necessary. We will then be able to directly 
observe possible quenching effects because in the comparison with [3] we can use exactly the 
same experimental input. While the agreement with [23] is not convincing at the moment, it does 
not seem worrying given the present errors. We emphasize that within about 100 MeV all these 
numbers agree, while a more precise statement needs additional work, in particular the continuum 
limit on our side, and on a longer term Nf > 2. 

3.3 B meson decay constant and Vub 

Let us now discuss /b. In HQET we consider the particular combination ^ 

<J>i = fBVm/2, (3.20) 
see eq. (3.3). We separate the static and the l/m^ contributions as 

log(rf Oi) = logirf^r) + [log(<I>i)](VH , (3.21) 

[l0g(<I.i)](V-) = [log(ZA)]('/-) + C<^in/^" + «spm/P" + (c2^-acr)/'^'''. (3.22) 

On the left panel of Fig. 7 we show that the GEVP method works as well as in the quenched 
case to extract the static matrix element. A good plateau is visible and our confidence in it is 
also based on the knowledge that corrections ai^e 0{e^^'^+^-"'°) when the computation is done as 
here (see [6, 15]). We have again applied two kinds of extrapolations of the static approximation 
^3/2jjjstat fj^g physical point [24]: 



r'J^^f' = b + anl (LO) , (3.23) 

(HM;t;PT) . (3.24) 



,3/2^stat 



3 1+3|2 2 



1 ^tni Inimz.) + c'riiz 



^We warn the reader that in this section we follow the notations introduced in [6], where the subscript 1 represents 
the ground state. In particular <I>i should not be confused with the finite volume meson mass <I>i introduced in section 2. 
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2 
1.9 
1.8 
dT 1.7 
'iH° 1.6 
1.5 



. N,. = 2, = 550 MeV 

. = 2, = 300 MeV 

« = 0, m = m 

» f ' l| s 



0.2 



0.4 



0.6 



0.8 



to [fm] 



1.5 



1.4 



■e 1.3 




b + C ITljj 

b'[l - 0.75(l+3g") m^" ln(m^0/(4n:f^)" + c" m^" 



1.2 



1.6 



Figure 7: On the left panel we show the plateaux of pf^' with A^=3 and t — to ^ 0.3 fm; the quenched 
data correspond to a lattice spacing quite close to the CLS setup at j3 = 5.3. On the right panel we show the 

3 /2 

chiral extrapolation of Oj'^'. 



Again we fixed g [7]. Tiie same fit form is used for Tq^^O^^^^. As sliown in tlie right panels of 
Fig. 7 and Fig. 8, wlietiier we do or do not include the chiral logarithm of HM;^PT, changes the 
value at the physical point by a small but noticeable amount. At the moment we take the average of 
the two extrapolations as the central value and include half of the difference as part of the systematic 
error. We obtain 

/bn?=2 = 178(16)(?).MeV, (3.25) 

where the first error includes the statistical uncertainty on matrix elements, the systematics coming 
from chiral extrapolation and the uncertainty on the physical scale of tq, while cut-off effects are 
not estimated yet. Let us recall that in the quenched approximation we have obtained Jb^^^q = 
234(18) MeV [26] for ro = 0.475 fm and the error covers for a 5% uncertainty on tq. 

At present there is a tension between the two ways of determining Vub from the exclusive B- 
decays B — )• Tilv and B — )• TV [27]. Both of them rely on the determination of the hadronic matrix 




Figure 8: Plateau of log(<I>i)(i/'^) - [log(ZA)]('/'") for N=4, t-to=4a, z^l5 and m,: ~ 550 MeV (left 
panel) and chiral extrapolation at the static order and at order 1/m (right panel). 
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Group 


Method 


Renormalisation 


Range in a 


ALPHA 
ETMC 

HPQCD 
FNAL/MILC 


HQET at 0(1 //Mb) 

twisted mass 
extrapolation in m\, 

NRQCD 

Fermilab action 


non-perturbative 
non-perturbative 

1-loop FT 
1-loop PT^ 


0.07 fm 
[0.065 -0.1] fm 

[0.09 -0.12] fm 
[0.09 -0.12] fm 



Table 4: Present methods in lattice computations of /b in unquenched simulations. 



elements from lattice QCD. In the Standard Model, the branching ratio for the latter decay is simply 
proportional to /glVubP with all other factors known. In Fig. 9, we have collected the values of /b 
obtained recently by several collaborations [25], from Nf = 2 (ALPHA, ETMC) and Nf = 2 + 1 
(FNAL/MILC, HPQCD) simulations. The heavy quark is discretised differently in each of the 
quoted computations: NRQCD (HPQCD), Fermilab action (FNAL/MILC), Twisted-Mass QCD 
(ETMC) and HQET (ALPHA). The systematic errors and difficulties in these different approaches 
are therefore rather different. For NRQCD a continuum limit does not exist and one studies a range 
of lattice spacings in the window between sizable discretisation en^ors and the divergent behaviour. 
In the relativistic twisted-mass approach one extrapolates in the heavy quark mass by assuming a 
HQET inspired form, and a continuum extrapolation has been performed. Some computations rely 
on perturbatively determined parameters. The FNAL/MILC and HPQCD computations use 1-loop 
renormalization^ for the currents and part of the parameters in the Lagrangian. In our approach, all 
renormalization is done non-perturbatively and the theory has a continuum limit, but so far only a 
single lattice spacing of a = 0.07 fm is available in large volume. There is an intrinsic truncation 
error of 0(A^/m^) which is present (more or less explicitly) in all approaches.^ We have sketched 
in Table 4 some of these issues. Despite the rather different characteristics, the results are very 
similar. The range of /b in Fig. 9 does not reach large enough values to remove the tension between 
B — > TV and B — nlv branching ratios within the Standard Model, and from Fig. 9 it seems very 
unlikely that the lattice determinations of /b are the reason for the tension. However, interpreting 
this as a hint for physics beyond the Standard Model, one has to keep in mind that the B — )• nlv 
form factors are less well studied than /b and that the experimental determination of the branching 
ratio B — )• TV is rather difficult. 

4. Conclusions 

We have reported on the status of the project undertaken by the ALPHA Collaboration to 
extract relevant B physics quantities from Nf = 2 lattice simulations in the framework of HQET 
expanded at 0(1 /m\,). The non-perturbative matching of HQET with QCD, through simulations 
performed in a small volume Li ~ 0.5 fm, is almost done. The measurement of HQET energies and 
matrix elements, using the GEVP approach, has started recently on ensembles produced by CLS. 
The analysis of quantities like the Z;-quark mass and the B decay constant is on the way. The first 

■'FNAL/MILC splits the renormalization factor into a non-perturbatively computed part and a rest estimated by 
1-loop perturbation theory. 

^Generically the truncation error is A^/rn^, but for the extrapolation in m significantly smaller masses enter. 
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Figure 9: Collection of recent lattice computations of /b . 

results are promising and once we have controlled cut-off effects by simulation at several lattice 
spacings, we will also determine hadronic parameters of the B — B mixing or B — )• tt semileptonic 
form factors as well as more details of the spectrum of hadrons with a b-flavor. 
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